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ABSTRACT:  Distribution  theory  for  likelihood  ratio  statistics  for 
the  comparison  of  several  treatments  with  a  control  is  dis¬ 
cussed.  These  test  statistics  account  for  prior  information 
that  the  treatments  are  at  least  as  effective  as  the  control. 


It  is  assumed  that  the  sample  sizes  on  the  treatments  are 
(approximately)  equal  and  the  sample  size  on  the  control  is  at 
least  as  large.  Normal  means  are  compared  under  the  assump¬ 
tion  of  a  common  variance,  either  known  or  unknown.  The  anal¬ 
ogous  problem  for  proportions  is  also  considered. 


Key  words  and  phrases:  Order  restricted  inference,  com¬ 
parison  with  a  control,  chi-bar-squared  distribution,  E-bar- 
squared  distribution. 
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1.  INTRODUCTION  AND  SUMMARY 

We  consider  an  experimental  situation  where  one  wishes  to 
compare  several  treatments  with  a  control  or  standard.  For 
example,  in  a  drug  study,  several  drugs  may  be  compared  with  a 
zero  dose  control.  In  the  absence  of  prior  information  about 
the  expected  responses  to  these  treatments  one  could  use  the 
results  of  Dunnett  (19551196^).  On  the  other  hand,  we  consider 
a  situation  where  the  investigator,  because  of  known  proper¬ 
ties  of  the  treatments,  wishes  to  carry  out  a  test  making  use 
of  prior  information  that  all  of  the  treatment  means  are  at 
least  as  large  as  the  control  mean.  (The  case  in  which  all  of 
the  treatment  means  are  no  larger  than  the  control  mean  is 
included  by  changing  the  signs  of  all  the  means.)  Assuming 

that  the  observations  are  normally  distributed  with  common 
2 

variance  o  ,  let  denote  the  hypothesis  that  4q  ■&  4^ ; 

i  =  l,2,**«,k  where  k  is  the  number  of  treatments,  4q  de¬ 
notes  the  control  mean  and  4^;  i  =  l,2,***,k  denote  the  treat¬ 
ment  means.  The  hypothesis,  H^,  is  a  special  type  of  order 
restriction  as  discussed  in  Barlow,  Bartholomew,  Bremner  and 
Brunk  (1972).  The  problem  of  testing  homogeneity,  HQ :  4q  =  4^ ; 
i  =l,2,***,k  when  the  alternative  is  restricted  by  a  partial 
order  is  discussed  in  Barlow  et  al.  (1972)  and  the  problem  of 
testing  a  partial  order  as  a  null  hypothesis  is  discussed  in 
Robertson  and  Wegman  (1978).  The  appropriate  null  hypothesis 
distributions  for  the  likelihood  ratio  test,  in  each  of  these 


problems,  are  mixtures  of  standard  distributions.  The  specific 

form  of  the  distribution  depends  upon  whether  HQ  or  the  order 
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restriction  is  the  null  hypothesis  and  upon  whether  or  not  o 
is  assumed  to  be  known.  However,  the  mixing  coefficients  in 
each  case  are  the  probabilities  that  the  maximum  likelihood 
estimates,  subject  to  the  restriction,  have  a  specified  number 
of  distinct  values.  These  probabilities  depend  upon  the  order 
restriction  and  upon  the  sample  sizes  and  they  can  be  difficult, 
if  not  impossible,  to  compute. 

The  partial  order,  H^,  considered  in  this  paper,  is 
termed  a  tree  in  Barlow  et  al.  (1972).  The  mixing  coefficients 
can  be  obtained  by  numerically  integrating  their  (3-38)  and 
using  their  recursive  relation,  (3.23),  but  this  approach  is 
quite  complicated  for  even  moderate  values  of  k.  For  the  case 
where  the  sample  sizes  are  equal  and  the  total  number  of  means 
does  not  exceed  12,  the  mixing  coefficients  are  given  in  Table 
A. 6  in  Barlow  et  al.  It  is  quite  common  to  have  significantly 
more  observations  on  the  control  than  on  the  treatments  and  we 
borrow  an  idea  from  Chase  (197*0  to  find  a  simple  approximation 
in  this  important  case. 

Williams  (1971,1972)  and  Chase  (197*0  considered  the  case 
in  which  the  researcher  Is  willing  to  assume  not  only  that  the 
treatment  means  are  greater  than  the  control,  but  that  the  order 
ing  among  the  treatment  means  is  completely  known.  Thus  the 
order  restriction  on  the  k+1  means  is  total.  One  application 
they  had  in  mind  was  the  comparison  of  increasing  dosage  levels 


of  a  drug  with  zero  dose  control.  We  follow  Chase's  approach 
and  obtain  approximate  critical  values  for  the  tree  ordering 
with  an  increased  sample  size  on  the  control.  The  limiting 
values  of  the  mixing  constants  are  found  as  the  sample  size  on 
the  control  becomes  infinite.  The  critical  values  based  upon 
these  limiting  constants  and  the  critical  values  based  upon 
the  equal-weights  mixing  constants  are  obtained.  The  approxi¬ 
mate  critical  values  are  constructed  by  interpolating  between 
these  two  values.  P-values  are  computed  by  interpolating  be¬ 
tween  the  P-values  obtained  from  the  equal  weights  and  the 
ones  obtained  from  the  limiting  values. 

Siskind  (1976)  and  Grove  (1980)  observed  that,  for  total 
orders,  the  mixing  constants  are  fairly  robust  to  moderate  vari¬ 
ation  in  sample  sizes.  The  same  kind  of  robustness  holds  for 
the  tree  partial  order  and  so  the  results  given  in  this  paper 
provide  reasonable  approximations  even  when  the  sample  sizes  on 
the  treatments  are  not  exactly  the  same. 


2 .  ONE-SIDED  TESTS 

In  order  to  be  specific,  consider  testing  the  hypothesis, 
Hq,  against  the  alternative  H^-Hq  (i.e.,  but  not  HQ) 

where  HQ  and  ^  are  specified  in  the  previous  section.  The 
data  is  obtained  by  choosing  k+1  independent  random  samples; 
applying  the  control  to  one  (of  size  ng)  and  a  different 
treatment  to  each  of  the  remaining  k  samples  (of  sizes  n^, 
n2*’  *,nk^*  Assume  that  the  resulting  sample  means. 


X.^  jO2/!!^);  1  *  0 ,1 ,2 ,  •  •  •  ,k. 


Accession  ? or 
KT1S  G?4l 
PTIC  TAP 

Ju.:1  c :: 


We  first  consider  the  case  in  which  a  is  known.  The 


likelihood  ratio  test  (LRT)  statistic  is 


Av::i  1 


T01 


•  C  ni^y 


where  fl  -  ^.0n1xi^,0n1  and  (Pq.Uj, •  •  •  ,1^) 


ft 


are  the  maximum 


likelihood  estimates  (MLEs)  under  the  restriction  H^.  Barlow 
et  al.  (1972)  discuss  the  computation  of  the  restricted  MLEs 
for  an  arbitrary  partial  order.  One  algorithm  that  can  be  ap¬ 
plied  for  any  partial  order  is  the  minimum  lower  sets  algorithm 
and  this  algorithm  takes  a  very  simple  form  for  our  tree  order¬ 
ing.  Of  course,  if  XQ  s  for  i  =  l,2,-**,k  then  4^  = 

for  i  ■  0,1,2, ••• ,k.  Otherwise,  arrange  in  increasing  order 
the  treatment  sample  means  (the  control  mean  is  not  to  be 
included)  and  denote  them  by  s  ^(2)  s  ”  ‘  &  ^(k)’  Next> 

find  the  smallest  positive  integer  j  for  which 

Aj  =  (nQX0  +^1=1n(1)X(1)  )/(nQ  +I^n(1))  <  x(j+i)  (1) 


where  the  symbol  n(-jj  is  used  to  denote  the  sample  size  assoc¬ 
iated  with  x(i)«  Such  an  integer  will  exist  unless 

s  ^(k)  and  ln  this  case  set  j  =  k.  Now,  4q  =  Aj  and 

the  restricted  MLE  for  the  treatment  mean  is  either  Aj  or  the 
mean 

treatment  sample Adepending  on  whether  the  treatment  sample  mean 
is  included  in  ^n(i)^(i)  or  not.  We  illustrate  the  algo¬ 


rithm  with  an  example. 
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Example .  Suppose  k  =  3 ,  nQ  =  25,  n-j_  =  n2  =  n^  =  10,  XQ  =  12.2, 
Xx  =  13.1,  X2  =  10.8  and  X^  =  11.9.  Since  the  sample  means  do 
not  satisfy  the  restrictions  in  we  need  to  apply  the  algo¬ 
rithm  to  obtain  the  restricted  MLEs .  Clearly  =  10*8, 

X^2j  =  11.9,  ^(3)  =  13.1  and  A.^  =  11.8.  Because  A ^  <  X(2)» 
j  =  1  and  =  jl2  =  11.8,  u1  =  13.1  and  =  11  *9  • 

Returning  to  the  test  being  considered  and  appealing  to 
Theorem  3.1  of  Barlow  et  al.  (1972),  we  see  that 

BCTqiXO2  2c]  ■  o<  /.k+DPCx^ic]  (2) 

2 

where  xv  denotes  a  standard  chi-squared  variable  with  v 

2 

degrees  freedom  (xQ  55  0)  and  P(-2,k+l)  is  the  probability, 
under  HQ ,  that  u  *  (4Q  •  •  •  ,4^)  has  exactly  l  distinct 

values.  The  P(i,k+1)  depend  on  the  values  of  no ,ni > * " ’ ,nk ’ 

In  general,  for  k  s  3  the  P(/,k+l)  can  be  obtained  by  the 
explicit  formulas  discussed  on  page  146  of  Barlow  et  al.  (1972) 
and  for  k  s  4  they  can  be  obtained  by  their  recursive  relation 
(3.23)  and  repeated  numerical  integration  of  their  (3.38).  Even 
for  moderate  k  this  may  require  several  numerical  integrations. 
If  nQ  *  n.^  ■  n2  *  •  •  •  =  n^  the  values  of  P(/,k+l)  are  given 
in  Table  A. 6  of  Barlow  et  al.  (1972)  for  k  £  11.  If  the  sample 
sizes  are  not  all  equal  but  do  not  vary  dramatically  (say  the 
ratio  of  the  largest  to  the  smallest  is  *  2)  then  the  values 
in  Table  A. 6  provide  a  fairly  reasonable  approximation. 

We  are  primarily  Inter°sted  in  the  case  where  the  sample 


size  for  the  control  is  significantly  larger  than  the  sample 
sizes  for  the  treatments.  In  the  special  case  in  which 
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n.^  *  n2  =  •  •  •  =  n^  and  n^/n^  is  a  positive  integer,  the  tables 
in  Ruben  (195*0  may  be  used  to  recursively  generate  the  P(i,k+1) 
(see  the  discussion  on  page  lM6  of  Barlow  et  al.  (1972)).  Even 
in  this  special  case  their  computation  may  be  very  tedious  so 
that  good  approximations  are  of  interest. 

As  was  noted  in  the  introduction.  Chase  (197*0  found  a 
good  approximation  for  the  case  where  nQ  is  significantly 
larger  than  n.^ ;  i  =l,2,***,k  but  when  specified  that 

s  ^  s  ...  s  He  assumed  that  n^  =  n2  =  •  .  •  =  n^  =  n  and 

obtained  the  limiting  critical  values  as  w  =  nQ/n  — ►  ®.  An 
Interpolation  between  the  equal  weights  critical  values  and 
these  limiting  values  worked  very  well  for  l  s.  w  <  ®.  We 
employ  an  analogous  approach  for  the  tree  ordering. 

The  P(2,k+1)  are  computed  under  HQ  and  so  they  depend 
only  on  the  variances  of  the  X ^ .  Because  the  weights  used  in 
pooling  the  sample  means  to  obtain  4  are  the  reciprocals  of 
these  variances,  it  is  common  to  let  P(/,k+l;  wQ ,w^ ,  • • • ,w^) 
denote  the  probability  of  exactly  l  distinct  values  in  il 
when  V(Xi)  88  w"1  for  i  *0,l,***,k. 

Theorem  1.  If  0  <  w.^  <  •  for  i  «l,2,***,k  then 
P(/,k+l;»,w1,w2,  •  •  •  ,wk)  -  lim  P(  i,k+l ;wQ  ,w1,w2  ,  •  •  •  ,wR)  «  2_k 

Wq  '  ' 

Proof:  Intuitively,  as  wg  — *  ®,  XQ  becomes  degenerate  at 


(3) 
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4q •  In  the  pooling  process,  any  time  one  of  the  treatment 
means  is  amalgamated  with  XQ  we  place  infinite  weight  on  XQ 
so  that  the  pooled  value  is  equal  to  Xq  .  Thus  P ( i,k+l ; “jW^ , 

•  •  •  ,wk)  is  equal  to  the  probability  that  exactly  1,-1  of  the 
treatment  means  are  not  amalgamated  with  Xq  (i.e.,  exactly 
i-1  of  the  treatment  means  exceed  the  common  value  of 
H0»U1J  *  '  '  »4k>  *  Thus>  P(i,k+l;®,w1,- •  •  ,wk)  = 

A  more  rigorous  proof  is  obtained  by  induction.  Obviously, 
for  k  *  1,  P(l,2;w0,w1)  =  P(2,2;wq,w;l)  =  1/2  for  all 

0  <  WqjW^  <®.  Assume  the  result  is  valid  for  k  =  m-1  and 
consider  k  =  m.  We  use  the  representation  for  the  P(i,k+1) 
given  in  the  recursive  relation  (3.23)  in  Barlow  et  al.  (1972) 
(which  is  valid  for  any  partial  order).  First  some  notation  is 
needed.  Let  k+1  be  the  collection  of  all  partitions  of 

0=  (o  ,  1 , • • •  ,k  ]  into  nonempty  sets  ,B2 , • • • ^  with  the  B^ 
sets  on  which  u  may  be  constant.  For  t  £  k,  such  decompo¬ 
sitions  are,  with  probability  one,  of  the  form  B^  =  £ 0 3  U  A 
where  A  c  {l,2,***,k3  with  card(A)  =  k-^+1  and  B2,**',Bj£ 
are  singletons  with  U  B.  =  {l,2,---,k]  -A.  Also  set 

J  ”  J 

WBi  =  Ej€B  WJ  and  for  B±  =  £j1<J2<**-<Jr>3  set 

w(B  )  =  (w  , w ,  ,-*-,w.  ).  For  1  <  l  s  m,  (3.23)  gives 
x  J  ^  J  2  Jp 

P(i,m+l:w0,w1, • • • ,wm)  = 

L,  P(/;/:WR  , • • • WR  )P( l,card(B1 ) :w(B, ) ) . 

Because  card(B1)  s  m  if  l  >  1,  we  take  the  limit  as  wQ  — ♦  ®, 


apply  the  Induction  hypothesis  and  note  that  card(^  m+^)  = 

{Jin)  -  (A)  t0  °staln 

P( /,m+l ; ®,w^ , • • • ,wm)  =  ^^^2  m  for  /=2,***,m. 

Next  we  note  that  by  ( 3 -  38)  of  Barlow  et  al.  (1972) 

P(m+l,m+l;w0,w1 ,  •  •  •  ,wm)  =  J*  ^-pz ~ j  >0(x)dx, 

where  *  (0)  is  the  c.d.f.  (p.d.f.)  of  a  standard  normal  dis¬ 
tribution  and  A^  =  w^/(Wg+w^)  for  i  *1,2,*  •  •  ,m.  As  Wq  ♦  “ 
A.^  — ►  0  and  so  P(m+1  ,m+l ;  “jW.^ ,  •  •  •  ,wm)  =  2  m  and 

P( 1 .m+1 ; “,w, , • • •  ,w  )  is  found  to  be  2~m  by  subtraction.  The 
'  *  ’1  m 

proof  is  completed. 

Table  1  contains  the  a  =  .1,.05  and  .01  critical  values 
for  the  statistic  Tq-^/cj  with  w  =  n^/n  =  1,®  and  k  =  2,3, 
••*,10.  The  w  -  1  values  were  obtained  by  using  (2)  with  the 
P(-£,k+l)  given  in  the  table  in  Barlow  et  al.  and  the  w  =  00 
values  are  based  on  the  P(«£,k+1)  given  in  Theorem  1.  For 
1  <  w  <  ®  we  recommend  interpolating  on  1  /*/w.  To  give  an 
indication  of  the  accuracy  of  this  approximation  the 
P(/,k+l;w,l,  •  •  •  ,1)  were  computed  for  k  =2,3,5  and  w=2,i4. 

The  true  significance  level  of  the  test  at  the  approximate 
a  =  .05  cutoff  value  (obtained  from  Table  1)  was  then  computed. 
The  largest  discrepancy  for  these  six  values  was  .0028. 


9 


Distribution  theory  for  normal  populations  provides  asymp¬ 
totic  theory  under  a  variety  of  assumptions  about  the  underly¬ 
ing  populations  (cf.  Theorem  4.5  in  Robertson  and  Wegman  and 
Robertson  (1978)).  For  example,  one  might  want  to  compare  the 
treatments  with  the  control  by  comparing  the  proportions  of 
individuals  who  exhibit  a  particular  response  to  the  stimuli. 

Let  p^  (pQ)  be  the  proportion  of  individuals  in  the  popula¬ 
tion  of  interest  who  will  exhibit  the  given  response  when  treat¬ 
ment  i  (the  control)  is  administered.  We  wish  to  test 

Ho:  P0  -  Pi  -  •••-  pk  vs-  Hi-Ho  wlfch  H1:  p0  ipi;  1  = 1’2>  ’  •  •  »k- 
If  ,p^  denotes  the  sample  proportion  from  a  sample  of  size  n^ 

and  if  the  samples  are  Independent  then  the  MLEs  of  the  p^ 
subject  to  can  be  obtained  from  the  algorithm  for  If 

the  are  replaced  by  p.^ .  Let  p ^ ;  i*0,l,***,k,  denote 

these  restricted  MLEs  (cf.  page  40  of  Barlow  et  al.  (1972)). 
Following  the  arguments  used  to  prove  Theorem  4.5  and  Corollary 
4.6  in  Robertson  and  Wegman  (1978),  we  see  that  the  LRT  statis¬ 
tic,  -2  In  A,  is  asymptotically  equivalent  to 

T01  =  [p(1-p)J-1  •  ^-oni^pi"p‘>2  as  nI  — *  °°  with  ni/no 

— — *  w1  €  (0,®)  for  i  =  0 ,1,2 ,  •  •  •  ,k,  where  p  =  E^^Wj^/E^  . 

Furthermore,  under  the  same  assumptions  on  the  sample  sizes, 

PCT^  ic]  —  Li_1  P(/,k+l)P[xMsc] 

where  the  P(/,k+l)  are  computed  with  respect  to  the  tree 
ordt  with  v  .ghts  wo ,wi» *  * ' ,wk’  Thus  T*^  is  a  test  statis- 
can  be  used  for  large  sample  sizes  and  if 


tlC  t  h  ct  o 
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=  n2  =  •  •  •  =  =  n  and  nQ/n  =  w  £  1  then  Table  1  can  be 

used  to  obtain  approximate  critical  values. 

We  now  return  to  testing  homogeneity  of  means  versus  the 

tree  alternative,  but  we  assume  the  common  variance  is  unknown. 

The  LRT  is  developed  in  Barlow  et  al.  (1972)  and  for  arbitrary 

sample  sizes,  the  test  statistic  is  a  monotone  function  of  the 

2 

ratio  of  estimates  of  a  under  the  two  alternatives.  Spe¬ 
cifically, 

S01  *  ILni(5i-'I)2<oI>i(xiJ-C)2 

where  the  (1^  are  the  restricted  MLEs  (computed  using  the  algo¬ 
rithm  described  earlier),  jl  is  the  weighted  sample  mean  (the 
sample  sizes  are  the  weights)  and  the  are  the  actual 

observations.  Theorem  3.2  of  Barlow  et  al.  gives  the  null 
hypothesis  distribution  of  Sq^  as  follows:  for  c  ^  0 

k+1 

P[S01  =  1  p(*£>k+1)p£B(i_ l)/2,(N-X)/2 

where  B  .  denotes  a  Beta  distribution  with  parameters  a  and 

cL  %  D 

b  (B_  .  =  0),  N  =  E  n,  and  the  P(/,k+l)  are  defined  as 

u,d  i=o  1 

before.  It  is  clear  that  nQ  — *  ®  implies  that  N  — *  <*>  which 
in  turn  implies  that  the  distribution  of  becomes  degener¬ 

ate  at  zero.  However,  Sq^  can  be  written  as 

n1(jli-u)2/(l^  Qni  (^-A)2  +  vS2)  where  vS2  is  distributed 

as  o2^  and  is  independent  of  the  sample  means  and  v  *  N-k-1. 


Following  Chase’s  work,  we  consider  the  distribution  of  the 

statistic  S*  =  Lk  n, (4  -4) 2/(Ik  n.  (X, -4) 2  +  Q)  where  Q 
01  i=o  1  1  i=0  1  1 

2  2 

is  fixed,  independent  of  the  sample  means  and  Q  ~  0  xv  (v  is 
now  a  fixed  positive  integer,  free  of  N).  The  same  arguments 
used  to  show  Theorem  3-2  in  Barlow  et  al.  give 

P[S01*c]  =  ^=1  p(A,k+Dp[B(je_i)/2j(k+v+i_je)/2  ic]. 

With  v  fixed,  n.^  =  n2  =  •  •  •  =  n^  and  w  =  n^/n^  — ♦  °°, 

P[S01SC]  -  2  ^i=l(i-l)P[B(i-l)/2,(k+v+l-i)/2  (4) 

Table  2  contains  the  a  =  .05  and  .01  critical  values  for 
this  limiting  distribution  of  S*^  with  w  =  nQ/n  =  1,“, 
k  =  2,3,  *  * ‘ ,10  and  v  = 2 , 4 , 6 , 8 ,10 ,12 ,14 , 16 , 18 , 20 , 30 ,50 , 100 . 

The  values  for  w  =1  were  obtained  by  using  the  P(-£,k+l)  in 
Appendix  A. 6  of  Barlow  et  al.  and  those  for  w  =  ®  were  obtained 
by  using  (4).  As  in  Chase's  work,  we  recommend  interpolating 
on  lA/w  for  1  <  w  <  “  and  on  N-1  =  (v+k+1)  1  for  inter¬ 
mediate  values  of  v. 

In  many  applications  of  the  above  procedures  the  researcher 
is  actually  interested  in  finding  evidence  in  favor  of 
rather  than  in  finding  evidence  against  Hq.  Testing  HQ 
against  H^-Hq  allows  one  to  control  the  probability  of 
falsely  confirming  H^.  If  one  is  more  interested  in  control¬ 
ling  the  probability  of  falsely  rejecting  he  also  might 

consider  the  test  procedure  developed  in  Robertson  and  Wegman 
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(1978).  They  considered  the  likelihood  ratio  test  where 

p 

is  the  null  hypothesis.  If  o  is  known,  the  test  statistic 
is  T12  =  -2  In  A  =  n^ilj-X^2.  Within  Hlt  HQ  is  least 

favorable  and  under  Hq  , 

P[T12  ^  PU  ,k+l)P[X^+1_je  ic]. 

Clearly,  if  ni>n2 *  * ’ ' »nk  are  fixed  and  n^  — *  ®  then 

P[T12  >d]  -  *0].  (5) 


The  equal  weights  critical  values  (taken  from  Robertson  and 

Wegman  (1978))  and  the  critical  values  for  nQ  *  00  (computed 

from  (5))  are  given  in  Table  3.  Again  we  recommend  interpola- 
-1/2 

tion  on  (n^n^  when  n^  =  =  •  •  •  =  n^  and  nQ  5s  n. 

2 

If  the  common  variance,  a  ,  is  unknown  then  the  likeli¬ 
hood  ratio  statistic,  for  testing  against  ~H^,  is  given 

by 


S 


12 


ni(Xi“Hi)2 

(Xi-Pi)2  +VS2 


and  the  appropriate  null  hypothesis  distribution  (computed  under 


P[S12  ic]  =  P(//,k+l)P[B(k+lije)/2j(N_k_1)/2  ic]. 

We  have  the  same  difficulty  we  had  in  testing  HQ  vs.  H^-HQ, 
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in  that  the  distribution  becomes  degenerate  as  nQ  (and  thus 

2 

N)  gets  large.  We  adopt  an  analogous  remedy  by  fixing  vS 


and  then  letting  n^/n^  — *  ®  with  n^  =  n^  =  •  •  •  =  n^.  Under 
these  conditions  we  find  that  the  appropriate  probability  is 


-k  yk+1/  k  \ 

vi=1U-iJ 


P[B(k+l ±Jt)/2,v/2 


Again  for  intermediate 


values  of  w  =  ng/n^  and  v  we  recommend  interpolation  in 
Table  4  on  1/Jv  and  on  N~^. 
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TABLE  1. 

Critical 

values  of 

TQ1/o2  for 

v  =  1,». 

a 

.1 

.05 

.01 

k 

w  =  1 

v  =  ® 

V  =  1 

w  =  00 

v  =  1 

v  =  ® 

2 

3.275 

2.953 

4.577 

4.231 

7.672 

7.283 

3 

4.696 

4.010 

6.171 

5.433 

9.561 

8.740 

4 

6.036 

4.955 

7.654 

6.500 

11.295 

10.020 

5 

7-333 

5.836 

9.075 

7.481 

12.939 

11.180 

6 

8.600 

6.672 

10.456 

8.4ll 

14.523 

12.275 

7 

9.848 

7.476 

11.821 

9.295 

16.061 

13.325 

8 

11.081 

8.257 

13.136 

10.153 

17.564 

14.300 

9 

12.301 

9.018 

14.446 

10.984 

19.039 

15.275 

10 

13.510 

9.764 

15.741 

11.800 

20 . 490 

16.200 
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TABLE  2  ( cont . ) 


V 

a 

wN 

k  2 

-  3 

4 

20 

0.05 

l 

.2029 

.2568 

.3013 

CO 

.1887 

.2288 

.2598 

0.01 

1 

.3164 

.3691 

.4121 

00 

.3027 

.3428 

.3721 

30 

0.05 

1 

,l4o6 

.1816 

.2166 

CO 

.1306 

.1611 

.1858 

0.01 

1 

.2246 

.2676 

.3027 

CO 

.2148 

.2471 

.2725 

50 

0.05 

1 

.0872 

.1144 

.1384 

CO 

.0808 

.1012 

.1183 

0.01 

1 

.l4l6 

.1719 

.1973 

CD 

.1353 

.1582 

.1768 

100 

0.05 

1 

.0447 

.0594 

.0727 

03 

.0413 

.0524 

.0620 

0.01 

1 

.0737 

.0906 

.1055 

CD 

.0703 

.0830 

.0938 

CD 

0.05 

1 

.0000 

.0000 

.0000 

00 

.0000 

.0000 

.0000 

0.01 

1 

.0000 

.0000 

.0000 

00 

.0000 

.0000 

.0000 

5 

6 

7 

8 

9 

_  10 

.3394 

.2849 

.3126 

.3062 

.4023 

.3245 

.4290 

.3406 

.4529 

.3547 

.4749 

.3672 

.4473 

.3955 

.4785 

.4l6o 

.5049 

.4316 

.5293 

.4463 

.5508 

.4590 

.5703 

.4697 

.2476 

.2065 

.2754 

.2249 

.3008 

.2410 

.3242 

.2554 

.3459 

.2686 

.3661 

.2805 

.3340 

.2930 

.3613 

.3115 

.3867 

.3271 

.4092 

.3418 

.4297 

.3545 

.4487 

.3652 

.1603 

.1332 

.1807 

.1466 

.1996 

.1588 

.2175 

.1702 

.2344 

.1805 

.2504 

.1904 

.2207 

.1924 

.2422 

.2010 

.2617 

.2191 

.2803 

.2314 

.2969 

.2422 

.3135 

.2524 

O  O 
— J  OO 

0 

4 er  M 

.0969 

.0784 

.1083 

.0857 

.1190 

.0926 

.1295 

.0991 

.1396 

.1055 

.1191 

.1035 

.1323 

.1123 

.1145 

.1206 

.1563 

.1282 

.1670 

.1353 

.1777 

.1421 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 
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TABLE  3. 

Critical 

values  of 

T12/o2  for 

V  =  1,® 

a 

.1 

.05 

.01 

k 

v  =  1 

w  =  00 

V  =  1 

v  =  00 

v  =  1 

w  =  * 

2 

2.5796 

2.953 

3.8232 

U.231 

6.8203 

7.283 

3 

3.1992 

u.010 

1*. 51*69 

5.1*33 

7.731*1* 

8.71*0 

1* 

3.6658 

U.955 

5.0830 

6.500 

8.1*082 

10.020 

5 

U.0U7U 

5.836 

5.5283 

7.1*81 

8.96U8 

11 . 180 

6 

1+.3630 

6.672 

5.8909 

8.1*11 

9.3926 

12.275 

7 

U.6U06 

7.U76 

6.2109 

9.295 

9.7969 

13.325 

8 

1*. 881*5 

8.257 

6.1*863 

10.153 

10.1250 

ll*.300 

9 

5.1013 

9.018 

6.7383 

10.981* 

10.1*297 

15.275 

10 

5.3013 

9. 761* 

6.9609 

11.800 

10.7311* 

16 . 200 
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TABLE  4  (cont.) 


V 

a 

\k  2 

3 

4 

5 

6 

7 

8 

9 

10 

' 

20 

0.05 

1 

.1760 

.2041 

.2241 

.2397 

.2524 

.2629 

.2720 

.2800 

.2871 

CD 

.1925 

.2373 

.2736 

.3045 

.3318 

.3562 

.3783 

.3985 

.4172 

0.01 

1 

.2920 

.3213 

.3418 

.3574 

.3701 

.3809 

.3896 

.3975 

.4043 

CD 

.3078 

.3527 

.3882 

.4l8o 

.4437 

.  4665 

.4869 

.5054 

.5223 

30 

0.05 

1 

.1206 

.1409 

.1558 

.1674 

.1769 

.1851 

.1919 

.1980 

.2036 

00 

.1324 

.1654 

.1928 

.2168 

.2383 

.2580 

.2762 

.2932 

.3091 

0.01 

1 

.2051 

.2275 

.2432 

.2559 

.2666 

.2754 

.2822 

.2891 

.2949 

ao 

.2169 

.2521 

.2807 

.3054 

.3273 

.3471 

.3651 

.3818 

.3973 

1 

50 

0.05 

1 

.0740 

.0870 

.0967 

.1042 

.1106 

.1160 

.1206 

.1248 

.1285 

1 

CO 

.0815 

.1029 

.1211 

.1374 

.1523 

.1661 

.1791 

.1914 

.2031 

j 

0.01 

1 

.1279 

.1436 

.1543 

.1631 

.1704 

.1768 

.1816 

.1865 

.1904 

( 

CO 

.1361 

.1601 

.1802 

.1979 

.2139 

.2287 

.2424 

.2553 

.2675 

I 

100 

0.05 

1 

.0376 

.0444 

.0496 

.0536 

.0571 

.0599 

.0625 

.0648 

.0699 

1 

00 

.0415 

.0529 

.0627 

.0716 

.0799 

.0878 

.0952 

.1024 

.109? 

0.01 

1 

.0662 

.0745 

.0806 

.0854 

.0894 

.0928 

.0959 

.0986 

.1011 

OD 

.0704 

.0836 

.0950 

.1051 

.1144 

.1232 

.1314 

.1393 

.1468 

00 

0.05 

1 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

00 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

0.01 

1 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

00 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 
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